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Abstract 

The one-electron density of states for the half-filled Hubbard model on a 
triangular lattice is studied as a function of both temperature and Hub- 
bard U using Quantum Monte Carlo. We find three regimes: (1) a strong- 
coupling Mott-Hubbard regime, characterized by a gap which persists even 
at high temperatures; (2) a weak-coupling paramagnetic regime, character- 
ized by the absence of a pseudogap at any finite temperature; and (3) an 
intermediate-coupling (spiral) spin-density-wave regime, characterized by a 
pseudogap which appears when U is increased beyond a critical (temperature 
dependent) value. The behavior of the y/3 x y/3 adlayer structures on fourth- 
group semiconductor surfaces is briefly commented upon in the light of the 
above discussion. 

Keywords Hubbard modeling. Semiconducting surfaces. Charge and 
spin density waves. Surface Mott insulators. 
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I. INTRODUCTION 

As is well-known, the frustration effects associated with the triangular lattice often lead 
to non-trivial ground-state degeneracies as in the antiferromagnetic (AF) spin 1/2 Ising 
model []I],H- The classical Heisemberg model on the two-dimensional (2D) triangular lattice 
with nearest-neighbor AF coupling and easy-axis exchange anisotropy is another example 
where frustration leads to a novel ground-state degeneracy. This system has attracted much 
attention especially since Anderson || suggested the possibility of a resonating-valence-bond 
ground state for the spin one-half case. Simply, the quantum liquid of radomly distributed 
spin-singlet pairs could be an efficient way to overcome the frustration of the Neel state in 
the triangular antiferromagnet. One more important source of interest in these lattices is 
the well-known diversity and richness of physical phenomena displayed by most transition 
metal compounds |§]. 

The experimental observation of low-temperature insulating phases in some v^3-adlayer 
structures on (111) Si and Ge surfaces has only fostered the interest in these 2D triangular 
lattices. Thus, whereas the -y/3-overlayers of Sn and Pb on Ge(lll) are both metallic at high 



1 



temperature, their corresponding low-temperature counterparts are either metallic, as in the 
case of Sn , or weakly insulating as in the case of Pb [[SHlOU . This latter system seems to 
go through some kind of reversible metallic to insulating transition whose precise nature is 
still controversial |TTHT3|. A charge-density wave (CDW) has been invoked in the case of Pb 
(but not in the case of Sn) as the driving force for the destabilization of the high-temperature 
phase, a conjecture not universally accepted at the time of writing. Related isoelectronic 
systems, on the other hand, like the \/3 x a/3 adlayer of Si on SiC(OOOl) [H, or of K on 
Si(lll):B [H| show a clear insulating behavior with a large gap and no phase transitions. 
These systems have been studied theoretically both within the local-density (LDA) []S|, [T~0|JT~T|] 



and the Hartree-Fock p~8| , [T9|] approximations. Quite recently |20| the LDA+U approach has 
been used to include strong on-site repulsions in SiC(OOOl). 

In this paper we carry out a model study of the triangular lattice in order to explore 
some general questions any realistic theory should comply with. For instance understanding 
the temperature behavior of the one-electron density of states (DOS) is essential to the 
development of a complete picture of these metal-insulator transitions. Hence we report the 
results of a Quantum Monte Carlo (QMC) simulation of the half-filled Hubbard model on 
such a triangular lattice in the grand canonical ensemble. The one-electron Green's function 
is studied as a function of both temperature and coupling constant (Hubbard U). As the 
temperature is lowered, a pseudogap develops in the one-electron DOS for intermediate 
values of U . This pseudogap is accompanied by two weak peaks in the spin structure factor 
which signal the formation of a spiral spin- density- wave (SDW). For lower U, no gap at all 
is found even for low temperatures, the system remaining always paramagnetic. For higher 
U, on the other hand, a well developed gap appears at any temperature, accompanied by 
a strong peak in the spin structure factor. The system is then brought into a state very 
similar to the ground state of the triangular antiferromagnet (the three-sublattice model). 
We emphasize that these are not distinct phases, but only different regimes with smooth 
transitions among them, as characterized by the behavior of the one-electron DOS. Since 
the presence or absence of a gap or pseudogap is of fundamental importance in determining 
the properties of a system, we believe that this type of characterization is useful and can be 
of help in understanding the electronic properties of the more complex adsorption systems 
referred to above. 

The paper is organized as follows: Section II recalls some of the basic properties of 
the triangular lattice as well as the Hubbard model, in order to fix the notation. Sec 
III summarizes very briefly the Quantum Monte Carlo algorithm in the grand canonical 
ensemble. The resulting one-electron DOS and spin structure factor are displayed and 
discussed in Sec IV and, finally, the paper closes with some concluding remarks in Sec V 



II. BASIC PROPERTIES OF THE TRIANGULAR LATTICE AND 

HAMILTONIAN 

In the -\/3 x a/3 R30° adlayer structures of Sn or Pb on Ge(lll) at one third coverage, 
each adsorbate sits on top of a triangle of Ge atoms. With just one unpaired electron per ad- 
sorbate, the overlayer is half-filled and, therefore, metallic in the absence of electron-electron 
interactions. This overlayer can in turn be described 3 lattice of adsorbate triangles 

(the three-sublattice model) with three unpaired electrons per triangle and, therefore, again 
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metallic in the absence of interaction. The corresponding surface Brillouin zones (SBZ) are 
the large and small hexagon, respectively, in Fig 1. The adlayer y/3 x \/3 lattice has just 
one band in the large zone, e° k , which folds onto three bands in the small zone. Fig 2 shows 
these three bands, e£, and ej~, unfolded in the extended zone scheme in order to see the 
nesting properties. They are given by (i=hopping strength) 

1 a/3 

e° k = 2t cos kx + At cos —k x cos ~^~^y (1) 

4 = ~^ £ k ± t ^> ^ sin k x -2 sin ~k x cos ( 2 ) 

It is easy to see that are just e° for k = (k x ± ^-,k y ). These bands cross at the points 
M' and K\ The wavevector if=(47r/3, 0) turns out to be a nesting vector with the band 
folding around the M'=0.5 K point (M'K — > MT, band e£). Likewise KM — ► TM' 
(band e k , ). Fig 3, finally, shows the resulting band along the small SBZ contour. Notice 
that the M'-fC' direction is obtained by folding the MK' portion of e°. 

This is of no consequence for the interaction-free system at half filling since the Fermi 
surface is not anywhere close enough to either the large or the small SBZ boundaries. When 
the interaction is turned on, however, the nesting symmetry may come into play, although 
weakly, at the points M and M\ closest to the Fermi surface. We shall see that, even at 
half-filling, this is indeed the case for the spin structure factor when U is large enough. 

In order to describe the interacting system, we adopt the Hubbard model, given by the 
standard Hamiltonian 

H = t Yl c ^ s ~ ^^2 n is + u Yl ~ 2^" iT ~ 2 V ^ 

<ij>s is i ^ ' 

where t is the hopping strength, U the on-site repulsion and fi the chemical potential. The 
single sums run over all the N x N adlayer atoms and the symbol Omeans summation 
over nearest neighbors (nn). As usual, cf s creates, while q s destroys, an electron of spin s at 
site i with occupation number = c^ s Cj S . We take £=0.055 eV so as to start with a narrow 
adlayer bandwidth (W=9t) of around 0.5 eV at U—0. U is varied to cover different regimes 
of the triangular lattice and \x is adjusted so as to have always half filling. Recall that, unlike 
the case of bipartite lattices, \i — U/2 does not necessarily correspond to half-filling since 
particle-hole symmetry does not hold in a triangular lattice. 

This Hubbard model is now simulated by the QMC approach in the grand canonical 
ensemble as initially developed by Blankenbecler et al and supplemented by a discrete 
lattice version of the Hubbard-Stratonovich transformation by Hirsch [J22[|. The whole ap- 



proach has been explained at length by Hirsch |23| and White et al |24[ and is briefly 
summarized in the following section. 



III. QUANTUM MONTE CARLO ALGORITHM 

In a grand canonical simulation |2"T|-F2~i| the imaginary time is discretized through the 
introduction of L time slices separated by an interval Ar such that (3=AtL. The partition 
function is then written as 
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Z = Tre 



-f3H 



Tre 



-AtLH 



Tr(e- ATK e- ATV ) L , 



(4) 



where 



C is C js 



}Z nis = zZ K 



*3 C is C js 



(5) 



<ij> 



<ijs> 



and 



(6) 



The last step in Eq(4) follows from the Trotter formula |25] which introduces systematic 
errors in the measured quantities of order At 2 . One should therefore like to take At as 
small as possible, although keeping the number of time slices not too large. This poses a 
serious limitation at low temperatures (large j3). 



A discrete Hubbard- Stratonovich (HS) transformation [2^] is now performed for each 
on-site interaction term within any time slice. 



(7) 



where a discrete, auxiliary Ising field an = ±1 (/ runs through the L time slices) has been 
introduced and coshAAr = e UAl ~/ 2 . Tr a means summation over the NL Ising "spins" an, 
which define the auxiliary Ising field through the lattice sites and time slices. Defining now 
N x N matrices V(l) such that 



we can write 



NL 

Z = ( V ArC/ / 4 ) Tr„TrlLjD;. 



where 



Df 



- AT T. ijs K H c ts c jsp- AT T.is Vu(i)a s ni. 



(8) 



(9) 



(10) 



where a s =±l for s=t i respectively. Notice carefully the difference between a, Ising spin, 
and s, fermion spin. 



The fermion degrees of freedom can be traced over [pl] , p3| to yield 

Z = Tr a U s det S [a] 
where S [a] , the fermion matrix, is given by 

s = I + B?Bf_ 1 ...B s l 

with 



(11) 



(12) 
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£>s _ e -ArK e -Ara s V(l) 



(13) 



and I being the unity matrix. Therefore IT S det S plays the role of a Botzmann weight for 
evaluating the average of any operator, < A >=TrAe~ l3H /Tre~@ H . The above steps must 
now be repeated for TrAe~^ H .Thus for A=Ci s Cj s , the resulting expression is 

Tr (c is c%e-P H ) = Tr a (Oj 1 ) .. II, det S [a] (14) 

and hence an equal-time single particle Green's function < Q s c+ > is obtained by averaging 
the corresponding matrix elements of the inverse of the fermion matrix S . Time dependent 



Green's functions can also be calculated and involve slightly more complicated matrices |23 
For any HS field configuration, Wick's theorem for operator products applies and one can 
therefore easily calculate time- dependent correlation functions in terms of single-particle 
Green's functions for any given Ising configuration. These imaginary-time quantities are 
then to be analytically continued to real time and frequencies via the maximum-entropy 
method [ 26 |27 which yields the corresponding spectral weight functions. 



The heat-bath algorithm is used to sample the HS field subject to the Boltzmann weight 
given above. This sampling is normally accomplished by single spin-flips an — >-<Ju, the 
corresponding weight-ratio being related to equal-time Green's function |2T| , |24|| . If the new 
configuration is accepted, the corresponding Green's functions is updated through simple 
operations |2"T| , |2"4| . Calculating the weight ratio requires N 2 operations, and therefore N 3 L 
operations must be done for a sweep through the whole Hs field. Normally two hundred 
warm-up sweeps and one thousand measurements separated by two sweeps were performed 
for each set of parameters. At low temperatures the algorithm becomes unstable due to the 
large number of time slices required in order to make At small enough. The product of B 
matrices in Eq(12) is ill-conditioned and one has to resort to matrix factorization techniques 



24 . 



IV. THE THREE COUPLING REGIMES 



According to the Mermin- Wagner theorem [28|, infinite-range magnetic order is forbidden 



in two dimensions at any T^O. This is so because the Goldstone modes strongly disorder 
the system giving rise to a spin-spin correlation length which decays with temperature as 
£(T) ~ exp(A/T), where A is a constant. Thus no phase transitions of magnetic origin can 
take place in an infinite system except, perhaps, at T=0. Other kinds of phase transitions 
are outside the scope of this theorem. Such is the case, e.g., of the (Mott) paramagnetic 
metal - paramagnetic insulator transition. Let us specialize to the case of the triangular 
lattice. 



A. The one-electron DOS 

The one-electron DOS is given by 

k 



where N is the number of lattice sites and A(k,uj), the spectral-weight function, is the 
imaginary part of the retarded one-electron Green's function. 

Mean-field studies at T=0 |29[ have shown that the half-filled triangular lattice is a para- 



magnetic metal in the weak-coupling regime, in contrast with the SDW insulating character 
of the square lattice for small U/t. No gap in the one-electron DOS is, therefore, expected 
at any temperature for an infinite triangular lattice. It has been shown, however, that size 
effects are very strong in this regime A gap in the one-electron DOS develops as soon 
as the spin-spin correlations extend over the whole system. Thus, for lattices of increasing 
size N x N, the system evolves from a situation where the correlation length £(T) > N, 
(with a gap) to one where £(T) < iV (without a gap). One should be careful when drawing 
conclusions about the existence of gaps from small lattices. 

Fig 4 shows the one-electron DOS of a half-filled 4x4 triangular lattice with periodic 
boundary conditions. We have taken U/t =5 (~ half the bandwidth, 9t), which is a weak 
to moderate value, and several values of /3, /3t=5, 10, 15 and 20. Even for fit as high as 
20, the system is far from having a fully developed gap. Since for the bigger lattices one 
expects weaker pseudogaps, it may be safely concluded that in the weak coupling regime a 
triangular lattice has no gaps at any temperature, in accordance with the Mermin- Wagner 
theorem. 

In the strong-coupling regime at T=0 the system is brought into a commensurate, three- 
sublattice, 120° twist SDW state (similar to the ground-state of the classical antiferromagnet) 
which is insulating and stable for increasing U. Quantum fluctuations about the classical 
antiferromagnetic solution lead to the essential qualitative physics of the Mott-Hubbard 
insulator at finite temperatures with a charge gap of order U in the spectral- weight function. 
Fig 5 bears the same information as Fig 4, but with U/t=20 which is deep inside the strong- 
coupling regime. Since size effects are very small in this regime JJT], it is fairly clear that a 
fully developed gap is present at any temperature. 

We thus see that, for a given temperature, the system evolves from a gapless situation 
at small U to a fully developed gap at large U. As U increases through the intermediate- 
coupling regime, one should find a critical value U C (T) for which the gap first appears. Fig 
6 displays, as Fig 4 and 5, the one-electron DOS for an intermediate value of U/t— 10 (~ 
the bandwidth). As the temperature is lowered from j3t= 5 down to f3t= 20, an incipient 
pseudogap gradually evolves into a fully developed gap. This value of U is clearly below 
the critical U for all 0t<2O, i.e., U c =10t for /3t= 20. The complementary view is given in 
Fig 7, which shows the one-electron DOS for (3t= 5 and U/t=5, 10, 15, and 20. We see the 
system evolving from a gapless regime to a pseudogap, a deep pseudogap and finally a fully 
developed gap. Thus U c —20t for (3t= 5. In this way one generates a temperature-dependent 
critical value of the coupling constant U C (T). 

B. The spin structure factor 

The spin structure factor, s(k), is given by the k- Fourier transform of the static spin-spin 
correlation function 

Sij =< (Ti z (T + )a jz (r) > r=0 (16) 
where r + =r + o + in the imaginary-time domain, and <Ji Z =n,^ — n^. 
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The peaks of s(k) and corresponding widths in /c-space give information, as is well- 
known, about the SDW's sustained by the system and the spin-spin correlation length. Fig 
8 shows s(k) for (5t— 20 and U/t—5, 10, and 20, representative values of the three regimes. 
The almost featureless shape for U/t—5 evolves into two sharp, although small, peaks close 
to M and M' in the intermediate regime and, finally into a large peak at M' in the strong- 
coupling regime. The system, correspondingly, evolves from a paramagnetic metal, through 
an incommensurate spiral SDW, into a Mott-Hubbard insulator 

V. CONCLUDING REMARKS 

The variation of the one-electron DOS with both temperature and coupling constant 
seems a useful tool for the purpose of identifying the different regimes of a given system. 
For the special case of the half-filled repulsive Hubbard model on a triangular lattice, we 
have identified an intermediate, temperature dependent coupling regime which interpolates 
smoothly between the weak-coupling (paramagnetic metal) and the strong-coupling (Mott- 
Hubbard insulator) regimes. As the temperature is lowered in this intermediate-coupling 
regime, the system evolves from metallic to insulating. 

We conclude with a comment on the \f?> x a/3 adlayer structures on group fourth semi- 
conductor surfaces. Although a close connection with the above model study is not claimed, 
these structures may constitute a physical realization of the three coupling regimes just 
described, Sn/Ge, Pb/Ge, and SiC being examples of the weak, intermediate and strong 
coupling regimes, respectively. Despite the added complexity due to electron-phonon in- 
teractions and atomic relaxation of both adsorbate and substrate atoms, the model study 
carried out here provides a general framework for the study of those systems. 

Acknowledgments. This work was supported by the DGICYT (Spain) Project 
N°PB98-0683. 
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VI. FIGURE CAPTIONS 



Fig 1. Large (outer exagon) and small (inner exagon) surface Brillouin zones 
(SBZ) of the triangular lattice. Shown are the especial points T, M', K, M, and 
K" 1 which delimit the contours Ti^Mrand TM'K'T used in the text. 

Fig 2. Band structure of the triangular lattice for U=0 in the three-sublattice 
model. The three bands, e\ (main band) and e^, are displayed along the TKMT 
contour of the large SBZ in order to show the band crossings and nesting sym- 
metry. 

Fig 3. Same as Fig 2, but along the TM'K'T contour of the small SBZ. 

Fig 4. One-electron density of states (DOS) of the triangular lattice for U/t—5 
(weak coupling) and decreasing temperature, (3t—5, 10, 15, and 20 

Fig 5. Same as Fig 4, but for U/t— 20 (strong coupling). 

Fig 6. Same as Figs 4 and 5, but for U/t— 10 (intermediate coupling). 

Fig 7. One-electron DOS of the triangular lattice for increasing U/t—5, 10, 15, 
and 20 at a fixed temperature (5t—5. 

Fig 8. Low-temperature (f3t=20) spin structure factor, s{k), of the triangu- 
lar lattice in the weak (U/t—5 ), intermediate (U/t— 10) and strong (U/t— 20 ) 
coupling regimes 
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